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A new class of solutions for Bose crystals with a simple cubic lattice consisting of N 
atoms is found. The wave function (WF) of the ground state takes the form \&o = 

e S™ +s b rj{ s i n (ki^Xj) sin (ki y yj) sin (ki z z 3 )}, where e Sb is the ground-state WF of a fluid, S$ = 



j 

q+k^O JV 



X>Sf ) (q)p- q + E S'^ 1) (q,k)p k p_ k _ q + . . ., p k = -±= J2 e lkr J , and k( = (ir/ai, n/ai, ir/ai) (oi 

q^O q.k^O V j = l 

is the lattice constant). The state with a single longitudinal acoustic phonon is described by the 



WF * k 



*l/o, where the 



P-k + X) ^(k, q)iV 1 Pqp_k-q + X Qi(g a; )p_k-q ;E + 7 permutations 

q 8as 

permutations give the terms with different signs of components of k. The structure of ^ k is such 
that the excitation corresponds, in fact, to the replacement of k; in some triple of sines from 
by k. Such a structure of and ^ k means that the crystal is created by sound: the ground state 
of a cubic crystal is formed by N identical three-dimensional standing waves similar to a longitu- 
dinal sound. It is also shown that the crystal in the ground state has a condensate of atoms with 
k = k; . The nonclassical inertia moment observed in crystals He 4 can be related to the synchronous 
tunneling of condensate atoms. 

PACS numbers: 61. 50. Ah, 67.80.-s, 67.80.bd, 67.80.de 



I. INTRODUCTION 

The science on crystals is developed for many years, 
and most properties of crystals are successfully explained. 
The current interest is focused on the regions of un- 
conventional and yet unguessed properties manifesting 
themselves, in particular, at extra-low temperatures (see 
surveys 1 '- 3 ). For the crystals with a charged lattice and 
the Bose crystals, such regions are, respectively, high- 
temperature superconductivity and supersolid phenom- 
ena. The splash-up of interest arose after the excellent 
experiments by E. Kim and M. Chan 4 ^, where a nonclas- 
sical inertia moment (NCIM) of a crystal He 4 was discov- 
ered. Later on, a number of new interesting properties 
joined by the term "supersolid" were found^— . There 
are almost no doubts that the "supersolid" phenomenon 
is related to the superfluidity of quantum crystals, which 
was predicted long ago^ 7 -. However, the physical nature 
of the superfluidity and the "supersolid" phenomenon is 
not clear yet, though a lot of models were developed^— . 

The results obtained below indicate that the basic 
property of crystals, i.e., the nature of crystalline order- 
ing, is not completely clear as well. 

Commonly accepted is the following structure of the 
WF of a Bose crystal^^: 

*o = e «>' IJ^ri-Ri), (1) 

i=i 

where N is the total number of atoms of a crystal, Yj and 
Rj are coordinates of atoms and sites of the lattice, the 
exponential function is the Bijl-Jastrow function taking 
correlations into account, and y(r) is usually written in 
the approximation of small oscillations: </?(r) = e~ a r I" 1 . 



By the modern ideas, a crystal is formed since its energy 
is less than that of a fluid. 

In what follows, we propose the basically new wave 
solution for the WF of a Bose crystal. In Ref. HH (cited 
below as I), it was shown that the states of a system of 
N interacting Bose particles positioned in a rectangular 
box L x x L y x L z in size include the states with the WF 

JV 

*o = Ae s » +Sb JJ{sin (ki x Xj) sin (k y yj) sin (h.Zj)}, 

3=1 

(2) 

where {k x ,k y ,k z ) = ki = (l x ir/L x , l y Tr/L y , l z n/L z ), 
l x ,ly,lz are integers, and the remaining designations 
are given in Sec. 2. In I, we analyzed only with 
ki = ki = (Tr/L x ,ir/L y ,it/L x ) (here, l x ,l y ,l z = 1), which 
describes the ground state of a gas and a fluid. In this 
case, the sines in @ form a standing half- wave that cov- 
ers the whole system and rests on the boundaries. But, 
l x , ly and l z may obviously take in the solution any other 
integer values except zero. It is natural to assume that 
if the half- wave is equal to the lattice constant, then WF 
(J2) describes a rectangular crystal lattice. Moreover, @ 
is one of the exact solutions of the Schrodinger equation 
with zero boundary conditions (BCs). In what follows, 
we will study this solution and the solution for a longitu- 
dinal acoustic phonon. We will show that the solutions 
agree with observable properties of crystals and predict 
a number of specific features, in particular, the conden- 
sates of phonons and atoms in a crystal. Solution @ 
testifies to the wave nature of a crystal. I did not deal 
with crystals earlier and have found the solutions acci- 
dentally, while studying the microstructure of a fluid. 

A short announcement of the results will be published 
separately 34 . 



2 



II. GROUND STATE OF A CRYSTAL 

Due to the presence of a product of sines in ©, 
the system can be partitioned into l x x l y x l z identical 
domains separated by plane surfaces, on which "So = 0. 
We may assume that if each domain contains one atom, 
and the size of domains is close to the equilibrium in- 
teratomic distance, then the system is stable. A crystal 
corresponds, obviously, to a system of domains with the 
size at which the energy of the system is minimum. If 
the number of domains is equal to the number of atoms, 
then l x = L x /R x , l y = L y /R yi l z = L z /R z , where R x , 
R y and R z are the periods of the lattice along the appro- 
priate axes. For a cubic crystal, they are equal to R. Let 
us study the properties of such system and make clear 
whether they correspond to the properties of crystals. 
We consider only a simple cubic (sc) lattice. 

Consider WF In it, Sb is given by the formula (see 

I) 



(tt) 



5 \ - «2(ki) 
6 = — 01 — P^iP-ki 



ki^O 
(tt) ki+k 2 #0 



+ E 



a 3 (ki,k 2 ) 



ki,k 2 #0 



3WN 



PkiPksP-ki-ka + • • ■ , (3) 



where 



/V 



3=1 

and the summation is carried on over the wave vectors 



L X Ly L z 



(5) 



Ux,jy:jz are integers), which are multiple to n/L. This 
is denoted by the symbol (tt) above the sums (under 
cyclic BCs, k are multiple to 2ir/L, and the solutions 

are strongly changed). The function S$ in ([2j takes the 
form of an infinite series (see I) 

q+ki#o (i), . 
cW oWc ^ 5 2 (<L k i) 

Lr\ „ l Lr\ \ I" 



q#0 



q,ki#o 



q+k!+k 2 #0 (1), , 

^ q, ki,k 2 ) 

2^ Jj Pk 1 Pk 2 P-k 1 -k 2 -q + --- (6) 



q,ki,k 2 #0 



where kj run values ([5]), and q take values 
2n (jf-, I 2 -, j^-j • In ([3]), the first sum is the Bijl-Jastrow 
function written in the variables p^. 

Let the faces of the crystal be ideally plane and paral- 
lel to atomic planes of the lattice. Every face creates a 
potential barrier for atoms of the crystal. We model this 
barrier for the face that is perpendicular to the X and 



has the coordinate x = by a step 
U w (x) 



U, 




x < 0, 
x > 0. 



(7) 



The potential of the face with the coordinate x = L x 
is U W (L X — x). Analogously, we can consider four other 
faces. For simplicity, we take U w = oo, which corre- 
sponds to zero BCs. 

In I, the product of assigning sines 

N 

* b s a c re (N) = l[{w(k, m x j )aii(k ly y j )8in(k,,z j )} (8) 

3=1 

has been factor out from the equations for and "Sk- 
in this case, we obtain the cotangent cot (irl x x/L x ), 
which cannot be expanded in a Fourier series, because 

J | cot (~kI x x/ L x )\dx = oo. To overcome this prob- 
o 

lem, we will use the following. It is easy to see that, 
under the change sinfca; — > |sinfcz|, the derivative 
(d\sinkx\/dx)(l/\smkx\) gives again cot (kx). There- 
fore, we replace *^ re © by the WF 

N 

* b s a c re (N) = I] {I sin {h xXj ) | x | sin (k ly Vj ) | x | sin (k lz zj) \ } 

3=1 

(9) 

and pass to 



/V 



(10) 



3 = 1 

x (\sm(ki y yj)\ + 5) x (| sin (hi,Zj)\ + 6)} , 



where 6 > 0, 5 — > 0. For such ^g" re , we obtain the func- 
tion 

f(x) = (df/dx)(l/f), f = \8m(k lm Xj)\ + 5 (11) 

instead of cot(fc/ x x). The former can be expanded in a 
Fourier series. For the singular point Xo, the series gives 
the arithmetic mean value of those at the points xq — 
and Xq + 0. In the final formulas, we transit to the limit 
5 — >• 0, which returns us to functions ([9]) and <JSj> - Bearing 
this fact in mind, we will immediately use in formulas the 
expansions at 6 — 0: 

cot (ki x x) = cot (irl x x/L x ) = Cl. (feK 9 * 3 (12) 
Here, q x — 2nj x /L x , j x runs all integers, and 
C lx (q x )=C lm (j x ) = 



i for j x — l Xl 2/^, 3l x , . . . 
i for j x = l Xl 2^, • 
for the rest j x . 



(13) 

The proof of formulas (JT2]) and (|T5|) is given in Appendix. 

We now have all the required in order to write the 
WFs of a Bose crystal. The WF of the ground state 
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is set by formula ([2]), where Sj and Sw are given by 
formulas ([3]) and ©• The equations for the functions 
a; and the ground-state energy of an atom (E ) follow 
from equations in I with the changes ki — > ki = k c , 

d(g) -> C,(g), Sj x) -> Sf : 



can be obtained analogously from the equations for S 
(see I): 

sP(q)eo(q) = 



(i) 



-4, 



2m 



(2tt 



(14) 



(15) 



q#0 



(27T) 



— E 2fc c , fe G x ( 9;c )Sf ) (-q x ) + (a; -> j/, z) 



-iVN2k cx q x a 2 (-(l x )Ci x (q x )Sq t q :i: + (x -> y, z) 

- £ {2(^ + q 1 q)4 1) (q,q 1 ) + 6^4 1) (q,q l! -q 1 ) 
qi#o 

^SPCqJS^q-qj.qx) 
JV^-qiqJS^Cqj^Cq-qx)} 

(2tt) 

E 2ik cx C lx (p x ) (q x -p x )S ( i\q-p x ) 



s* = 1LlL 

2 d + 1 AT 



(tt) 



(tt) 



2 d h 2 
(tt) 



1 1 ?i 2 £r 

™^(0)-9*HAf E E ^° 2 ( k )' 

(16) 
(17) 



2px4 1) (q-p 3 ,,-q) 



4%, q 1 )[e (q 1 ) + e (q + q 1 )] 



(22) 



+ 25} 1) (q)a 2 (-q 1 )qqi-g 2 5W(q)a 3 (q ) q 1 ) 
= Sq^VNikcxQ^qx) {2g lx a 2 (q 1 ) - g 2; a3(q 1 , qj} 

(27T) 

+ E ik cxCi x (p x ) 4(q x + q lx -p x )S!p(<l- Pj^qx) 

- e^iSPCq-Pa^qx^q-qx) + (a; y, z) 
1 00 f 2 

+ "Tt? E 1 "777^2 -qiq 2 )5 , 3 1) (q,qi -q 2 >q 2 ) 

4 (i) 

h — (qi + q 2 + q)"5 , 3 (q>qi.q2) 



2o0), 



a 2 (k)fc 2 -a^(k)fc 2 =A 2 (k) 

(tt) 



+ ^ E a3 ( k > *30C« a + k q) + ^ E « 2 ° 4 ^' k )' 



o 3 (k,q)« A 3 (k, q)-2 



q^O 



fl(k,q) 



eo(*:) +eo(g) + e (k + q)' 



(18) 



i?(k, q) = kqa 2 (k)a 2 (q) - k(k + q)a 2 (k)a 2 (k + q) 

- q(k + q)a 2 (q)a 2 (k + q), (19) 



1 (2?r) 

A 2 (k) = -= E ik cx C lx {-q x ) [4(q x + k x )S^(d x ,-k) 



" 4q 2 ( qi - q 2 + q)5< 1) (q 2 )^(q - q 2 , q x ) 
+ 6g 2 2 ^ 1} (q 2 )5«(q - q 2 , q^ -q - Ql ) 
+ 4(q 1 +q 2 ) 2 5( 1) (q 2 ,q 1 )^ 1) (q-q 2 ,q 1 +q 2 )}. (23) 

Equations (fT4")) -(f2"3" | are a complicated system of non- 
linear integral equations, whose solutions determine the 
properties of the ground state of the crystal with a rect- 
angular lattice. Let us analyze these equations. Of a 
paramount interest are the value of Eo and the distribu- 
tion of atoms in the crystal. 

It is seen from Eq. ([22"]) that the nonzero value of 

(q) is determined by the first "one-dimensional" term 
~ Cj*(fe)^q,q on the right-hand side. This yields the 

one-dimensional solutions 5^ (qJ^q.q,,,) &i (q tf )^q,q > 

and (q z )<5 q ,q, ■ However, the terms ~ Sip (q-JS'.^ (q _ 
qx , qx ) and some other ones on the right-hand side gen- 
erate also not one-dimensional solutions of the form 
Si\q x + %)Sq^+ci v and s[ 1} (q x + q y + qj of the 



(l)/ 



+ 6feS , | 1) (q c ,-q a .,k) + (x -> y, z) 

(2*0 

N 



1 ( ' r 

-E[% 2 + qk)5 1 1) (-q)5 2 1) (qk)+ 



q#0 



+ 4(k + q) 2 ^ 1) (q,k)4 1) (^q,-k) 
+ 6q 2 s[ 1) (-q)Si\q,k,--k.) 



(20) 



where eo(k) = fc 2 (l — 2a 2 (k)), d = 3 is the dimension of 
the system, and 



= / U 3 (r)e-^ T dr 



(21) 



is the Fourier transform of the interaction potential U3 (r) 
of two Bose particles. The equations for the functions 



same order (~ v N) as one-dimensional solutions. Since 
Cl x {jx) are nonzero only at the "resonance" points j x = 
±l x , ±2l Xl ±3Z X , . . ., the function 5^(q) is nonzero only 
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for the "resonance" wave vectors 

TPS/ \ r» ( Tx^X T z l z 

q res (r) = 27r — 

V J^x J^y J-'z 

where r = {r x ,r y ,r z ), r x ,r y ,r z = ±1,±2,±3, . 
regard for these relations, we can write the solution 

S[ i] (ci = q x + q y + qj = S^(q) {(J q , q ™ + <5 q , q j 



(24) 



With 



°q,q; es "f 
<W«> 



^q,q™ s +qj" + ^q™ +q; es +^q,qr: cs +qj 



(25) 



and, analogously, sj> 2 (q, q 1; ...q^j). As is seen from Eq. 
}, the values of ^'(qj.) depend significantly on the 



sums with S 2 \ Sfp and with the terms ~ S-p S± , and 

~ Si ' S 2 ^ ■ However, we restrict ourselves to the zero 
approximation 



W N2k cx Ci :c (q x )a 2 {-q x )5 cl;iiM r 
q x - 2q x a 2 (q x ) 



(26) 



1qi x a 2 (Qi) 



2a 2 (q x ) 
- fea 3 (q 1 ,q 2 . 



eo(qi) +eo(q K +qi) 



(27) 



The analogous relations are true for S^ 1 ' (q^), S^\q z ) 
and S<p(q y , q x ), (q z ,q 1 ). The omitted corrections 

can renormalize Si(q 3 .), by changing its value by sev- 
eral times. In this case, the higher corrections are 
damped by the decrease of a2(q) with increase in q. 
For the not one-dimensional solutions sf^q,,. + q y ) and 
•Si'^qx + % +q z ); even the zero approximation is a com- 
plicated sum (over q x ) of the terms ~ 5^' ) (q 1 )5'^(q— q x ) 

and ~ 5^(q 1 )52^(q — q 1 ,q 1 ). By our estimates, the 
not one-dimensional solutions are significantly less than 
the one-dimensional ones. Due to the complexity of the 
equations, we neglect the not one-dimensional solutions: 



sW(q)*S«(q){rf 



<W*}- (28) 



Let us study the distribution of atoms along the re- 
axis (coinciding with one of the axes of a crystal) for the 
lattice of He 4 atoms at R x = R y = R z = R = 3.6 A. 
According to I, the probability for an atom to be at a 
point x is approximately determined by the formulas 



■>jj(x) ~ sm(k cx x)e Sl{x \ 



si(x) 



(2tt) 

E 



S[ l) (q x ) 



N 



(29) 



(30) 



where j x = ±1, ±2, ±3, . . . and k cx = ir/R. Despite the 
approximate character of the formulas, we may expect 
that they give the general form of a probability distribu- 
tion in a cell. 

Since ^^(q^.) ^ only at the resonance points and 
^(-q^sftqj, we have 



si(x) 



E 



S[ l) (q x ) 



N 



2 cos (q x x) 



(31) 



where q x — r x 2nl x / L x — r x 2n/R. We will determine 
Si\q x ) from Eq. (|26p . by using the zero approximation 
for a 2 (q), 



a-2{q) 




nv^,{q)m 
8H 2 q 2 ' 



(32) 



which follows from (fTTj) if all sums are neglected. Note 
that we took the solution with the sign "minus" before 
the root (see I). 

We choose the interatomic interaction potential for 
atoms of the crystal, as in I: 



Ub r < a 
U bd a <r <b 
r > b. 



(33) 



Below, we use, unless otherwise indicated, a = 2 A, 
b = 4 A, Ubd = — 9K, which corresponds approximately 
to He 4 atoms. This potential is sufficiently crude, but 
it is qualitatively proper and has the analytical Fourier 
transform 



4-7T 



(34) 



[Ub — Ubd] [sin (ak) — ak cos (ak)] 

Ubd [sin (bk) - bk cos (bk)]. (35) 



By formulas P5 ]) -([55 ]I . we determine ip(x) (gSJ. It 
has a periodic shape, in correspondence with the domain 
structure. The distribution |^(x)| 2 in one of the domains 
is shown in Fig. [JJ The value of R is chosen like that for 
He II: R = 3.6 A, which is close to R of the crystalline 
phases of He 4 . As is seen from Fig. [TJ the distribution 
[■^(a;)! 2 at the height of the potential barrier Ub — 300 K 
is similar to the bare one (sin 2 (k cx x)), but is more flat- 
tened. At Ub > 1000 K, we see the appearance of two 
maxima located symmetrically relative to the center of a 
cell. They increase with Ub. So, the probability density 
is the highest not at the center of a cell, as is commonly 
accepted, but at these maxima. In the three-dimensional 
case, the maxima indicate the presence, inside a cell, of 
an orbit with cubic shape. The orbit depends on the val- 
ues of a and b: at a = 1 A, b = 3 A and Ub — 3000 K, the 
maxima disappear, and |^(a;)| 2 is similar to the curve of 
triangles in Fig. [TJ But such small a does not correspond 
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Fig. 1: Probability density |f/>(a:)| 2 lUSJl, (O versus the co- 
ordinate x of a cell of the sc crystal at the potential with 
a = 2 A, b = 4 A, U bd = -9K and the barriers t/ 6 = 300 K 
(triangles), t/ 6 = 3000 K (circles), and U b = 3000 K with re- 
gard for only the first term in sum (f3T|l (squares). Stars stand 
for the bare function sin 2 (k cx x). The z-axis coincides with 
the 2-axis of a cell of the crystal, values of x are given in A. 



Fig. 2: Function |i/>(:r)| 2 for a sc crystal of He 4 at high barriers 
U h . Potential ([33]) with a = 2 A, b = 4 A, U bd = -9 K and the 
barriers U b = 9000 K (crosses), U b = 10 5 K (squares), U b = 
10 6 K (circles) and U b = 10 7 K (triangles). Stars mark the 
curve for a crystal of krypton atoms with the same parameters 
except the mass (m = 2\m^) and potential (Ub ~ 3000 K, 
J7 M = -140K). 



to the He 4 -He 4 potential. As a increases by 1 A, there 
appears a clear orbit, like the curve of circles in Fig. [TJ 

The orbit size is approximately equal to a half of the 
cell size, because the main contribution to the maxima 
is given by the first term (with q x = 2k cx ) in sum (1311) . 
The terms with q x = Ak cx ,6k cx ,8k cx are small, and the 
peaks for q x = Ak cx arise only at U b > 10 6 K (see Fig. [5]). 
To calculate ^(a;)! 2 , it is sufficient to take two first terms 
in sum (J2TJ) (due to the fast decrease of a^iq), as q in- 
creases). 

In Fig. [3] we present the dependence of the orbit on the 
barrier height Ub- As is seen, the orbit becomes narrow 
at large U b . Moreover, at U b = 10 6 K and 10 7 K, the 
second orbit appears at a distance of ~ 0.5 A from the 
wall of a cell. However, the neglected higher correlative 
corrections become large at large Ub, which can lead to 
the widening of orbits. 

In Fig- El stars show ^(a;)! 2 for the lattice of krypton 
atoms. The orbit is narrow at the barrier Ub = 3000 K, 
whereas the orbit for helium atoms is wide at such U b - 
However, it is knownl>2£ that the ground state of the 
crystal of heavy inert elements is well described under 
the assumption of small oscillations of atoms near points 
of the lattice. Apparently, there is no contradiction in 
this case, since the function r 2 \ip(r)\ 2 at small oscillations 
is characterized by a spherical orbiti with approximately 
the same size. Nevertheless, the orbit in Fig.f2]has shape 
of the surface of a cube (with edge R/2), rather than a 
sphere, even for the function r 2 |?/;(r)| 2 . In other words, 
the motion of an atom is oscillatory only approximately. 
More exactly, this motion has a wave character. This 
fact is unusual and means that the resonance wave (the 



product of sines in ([5])) sets the lattice of a crystal and 
the motion of atoms in the cells. 

Can an orbit be discovered in experiments? The scat- 
tering of light or neutrons in a crystal will show ordinary 
Bragg- Wolf peaks first of all, since the time average of 
the positions of an atom is the center of a cell. But the 
orbit can be revealed in some specific features of the scat- 
tering. 

Let us estimate the ground-state energy Eq of a sc crys- 
tal, using formulas (ll~4l - (fT7|) . With regard for solution 

(J2SJ) for S[ l \q) with Sf^qJ ([H and solution $T5$ for 
Ci(j), we obtain 



E = E° + A 1 , A x = 



2m 



a2(q x ) ~ 3a 2 (q x ) 
(2a 2 (q x ) - l) 2 



(36) 

where r = 1, 2, 3, . . ., q x = 2nr/R, k 2 c = 3n 2 /R 2 , and E b 
(fT6| coincides with E of a Bose fluid. Consider a crystal 
of He 4 atoms with R = 3.6 A, h 2 k 2 /2m 13.85 K, and 
the potential (J33J) with a = 2 A, b = 4 A, Ubd = -9K. 
Using the zero approximation (|32p for a<i (q x ) and taking 
the barrier U b = 1000 K, we obtain At w -12 K, E w 
-430 K. For U b = 200 K, we have A x « 4 K, E -36 K, 
whereas ^ k 9K, £ - -IKfor [/ 6 = 81 K. The last 
value of Eq is close to the experimental onei. For He 
II, i?o corresponds to experimental data also for U b ~ 
100 K (see I). For the realistic value U b ~ 10 3 - 10 4 K, 
the experimental value Eq ~ —IK can be obtained with 
regard for correlative corrections. 

The most essential and unexpected is the conclusion 
that the lattice is created by a standing wave in the prob- 
ability field. As we see in Sec. 3, this wave is similar to 
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a sound one. In other words, the crystals have a wave satisfy the equations 
nature. 



III. STATE WITH A LONGITUDINAL 
ACOUSTIC PHONON 



Consider a crystal with sc lattice. An optical phonons 
are absent for it. We consider only longitudinal acoustic 
phonons. The WF of a crystal with a single phonon can 
be obtained from the solution for a one-phonon state of 
the Bose fluid (see I) by the following changes: ki — > k c , 

Ci(q) -> Ci(q), and sf 3 -> Sf\ We obtain 



* k (ri,...,rjv) = V>k*o, 



^Ak = + 7 permutations, 



(2tt) 

V k = V>k + bo(k)+ £ Q!(q,k)p_ q _k 

q#0,-k 

<32(q,qi,k) 

?TF PqxP-q-qx-k 

q,qi^o 
q+qi+q 2 +k#o 

E 

q.qi.q 2 ^° 



(37) 



(38) 



'JV 

Q3(q, qi, q 2 » k ) 

JV 



(39) 



" Pqi Pq 2 Z 3 - q- qi - q 2 - k 



(tt) 



v^ k = p-k+ 53 

k 2 ^0 : -k 
(tt) k 2 +k 3 +k^0 

E 

k 2 ,k 3 #0 



6 2 (k 2 ,k) 



JV 

& 3 (k 2 ,k 3 ,k) 



"Pk 2 /3-k 2 -k 



(40) 



TV 



Pk 2 Pk 3 P-k 2 -k 3 -k + 



where the permutation means ip^ with the different sign 
of one or several components of the vector k, the vector 
q is quantized like 2-irj/L, and the vectors kf, q-, and k 
are quantized like ftj/L. 



Solution (|38p describes a three-dimensional standing 
wave decaying into eight counter traveling waves. 

The energy of a phonon E(k) and the functions bj, Qj 



e(k) = eo(fc) - - b 2 (k 2 , k)2k 2 (k 2 + k) 



k 2 #0 



(tt) 



i 53 6fc 2 2 6 3 (k 2 ,-k 2 ,k) 



2 



k 2 ^0 

(27T) 

53<2 1 (q,k)(k + q) 

q#0 

q^ 1) (-q)+2(k + q)4 1) (-q,-k)' 

4 (27r) 

-7= ^QaC-* q, (q) 

ViV q^O 



(41) 



^ ] 2k cx (k x 

+ ^2 Ak cx q x iCi x (q x )Q 2 {-q x ,q x ,'k) + (x -> y,z), 

Qi(q,k) [e(fc)-e (k + q)] = 
: 2ik cx Ci m (q x )5q^ [-k x + 2q x b 2 (q, k)] 

2S'i (1) (q)r , „ 2, / , m 4fc 2 
1 VH Mqk- 2a 2 &9fq,k)l - - ; 

/JV 

1 

JV 



[qk - 2g 2 6 2 (q, k)] - ^S< 1} (q, -q - k) 



E {2qi(qi+q + k) 
qi#o 

x [Q 2 (q, q 1; k) + Qi(q + q 1; k) VnS? {-qj] 

+ 6<2 2 Q 3 (q,qi,-qi,k) 



WaCq-qi^kKJVSHqi) 
4( qi + q + k) 2 Qx(q + qi ,k)A/jv4 13 (- qi , -q - k)} 

(27T) 

53 2k cx iCi x {p x ) [{-k x - q x +p x )Qi(q-p x ,k) 
2p x Q 2 (q- P s! p x! k)] + (x -^y,z), 



(42) 



6 (k)e(fc) = -2fc 2 sf (-k)<5 k)k e 
2k cx k x iCi x (-k x )VNS k ,ki 
2 00 

-t= E < J 2 '?2(-k, q, k)(5k,k<= 



(2tt) 

53 2(k + q) 2 g 1 (q,k)^ 1) (-k~q)^ 

q^O ,-k 

(2tt) 

+ 53 2 ^ fc ca; feiC'^('Jx)Qi(-k- q a .,k)(5k,k<= 

+ (x-*y,z), (43) 
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6 2 (k 2 , k) [e(k) - e (k 2 ) - e (k + k 2 )] = 
kk 2 a 2 (k 2 ) - k(k + k 2 )a 2 (k + k 2 ) - fc 2 a 3 (k,k 2 ) 



(7T) 

E 

k 3 #0 



3k 3 (k 2 +k 3 + k) 



N 



6 3 (k 2 ,k 3 ,k) 



2(k 3 -k 2 )5 2 1) (-k 3 ,k 2 



Q 1 (k 3 ,k)(k 3 +k 

3(k 3 + k)5< 1) (-k 3 ,k 2 ,-k 2 -k 
2g 2 (-k 3 ,k 2 + k 3 ,k)(k 2 +k 3 ) 
k 3 ^ 1) (k 3 ) + 2(k 2 +k 3 )5^ 1) (k 3 ,k 2 ) 

3fc 3 2 ^ 1} (k 3 )g 3 (-k 3 , k 2 , -k 2 - k, k)} 



(2tt) 

E 2fc - 



dQJPx) [3p x Q 3 (-p a ,,k 2 , 

2(fc 2 a: +Px)Q2(-P x ,k 2 +p x ,k)] , 

higher corrections + (a; — > y, z), 



k, k) 



(44) 



where e(fc) = 2rnE(k)/h 2 , k e is the wave vector with 
all even components (i.e., they are multiple to 2n/L), 
q x = q x i x , and p x = p x i x . By (x y,z), we denote 
the same terms as one with the separated ir-component, 
but with the changes x — > y and x — > z. All terms 
of Eqs. pT j) - ([4"4")) . except for the first one, contain a 
product of two wave vectors (for example, q 2 , k cx k x , or 
(k 3 + k)(k 3 — k 2 )). These wave vectors must be nonzero. 

Some of these equations contain the sums, where the 
first argument of the functions Qi or Sj has components 

multiple to ir/L. In this case, we set Qi = 0, Sj = 0, 
because, by the definition of these functions, the compo- 
nents of the first argument are multiple to 2tt/L. 

Equations (|4T |) - p4|) are written in the approximation 
of "two sums in the wave vector", at which the series 
contain the functions a 2 , a 3 , & 2 , 63, Sj< 3 , and Qz< 3 and 

do not include the corrections aj>4, &j>4, «Sj>4j an d Qi>4- 
It follows from Eq. (|42p that the solution for the func- 
tion Qi(q, k) has a "resonance" form analogous to (1231) . 
We restrict ourselves to the one-dimensional approxima- 
tion: 



Q«(q,k)~Q«(q,k){5 q , qS ea 



(45) 

In the zero approximation without regard for the sums 
in Eq. (|42|) . we have 



(46) 



Ql(q x ,k) = [2ik cx a js (q x )-2q x S?\c ix )/VN 

-k x + 2q x b 2 (q x ,k) 
e(fc) - e (k + qj 
i2k cx Ci x (q x ) 



-k x + 2q x b 2 (q x ,k) 

l-2a 2 (q x ) A e(k) - e (k + qj 

The analogous relations can be given for Qi(q a ,k), and 
Oi(q z ,k). 



Equations (I4TT) — (|44[) are very complicated. By setting 
Qi = in them, we obtain the equations for a Bosc fluid, 
namely, for 6 2 and the dispersion curve E(k). Thus, the 
equations for a crystal are the equations for a fluid plus 
some additional anisotropic corrections. 

Let us consider dispersion curves. First, the value of 60 
does not influence E(k). In the simplest approximation 
with regard for the anisotropy, the dispersion curves are 
determined by the formula 

2 (27r) 

e(k) = e (fc) - -= Qi(q, k)(k + q)qSf(-q) 



(2*-) 



- ^2 2k cx{k x + q x )iCi l (-q x )Q 1 {ci x ,]<i) 

+ (x->y,z). (47) 

With regard for solutions (|T3" ]> , ([2"5 ]) . ([25 ]) . and (|g5 j) 

in the approximation 6 2 = 0, we obtain 

e(k) = £ „(fc) + e-(k) + e-(k) + ef (k), (48) 



^ r (k) 



^ (1 



k x (k x + (Ja; 



9x=gJ 



(l-2a 2 ( fe )) 2 e(k)-e (k + qj ! 



(49) 

where the sum is taken over the resonance values of 
q res _ ±2k CXl ±4k cx , ±6k cx , .... For a sc crystal, 

E(k , = ^ + M, ( 5 0> 

where j=l, 2, and 3 for directions (1,0,0), (1,1,0), and 
(1,1,1). Since the values of e£ r (k) and eo(k) are of the 
same order, the dispersion curves are different for differ- 
ent directions. In other words, the spectrum of longitu- 
dinal acoustic phonons is anisotropic, what was observed 
in experiments. 

Due to a large value of e^T, the method of iterations 
does not converge, and the more exact numerical meth- 
ods require much time for the analysis. To demonstrate 
the influence of the correction e x r , we decrease it by two 
orders of magnitude so that the method of iterations can 
be applied. In Fig. 3, we present the dispersion curves 
calculated by formula (|50l) . The curves are different for 
different directions and are similar to the observed ones 2 - 
for a bcc crystal of He 4 . 

The qualitative behavior of the dispersion curves can 
be understood, by using the simple formula 



E(k) 



H 2 e (k) H 2 k 2 



2m 



2m 



(l-2aa(k)). 



(51) 



With regard for approximation 
mula becomes 



for a 2 (k), this for- 



E(k) 



fh 2 k 2 



V 2m 



nv^{k) h 2 k 2 



(52) 
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Fig. 3: Dispersion curves E(k) (I50|l with e% r decreased by two 
orders of magnitude for directions (1,0,0) (triangles), (1,1,0) 
(circles), and (1,1,1) (squares) for a sc crystal of He 4 atoms 
with R = 3.6 A and potential J33]) with a = 2 A, b = 4 A, 
Ubd ~ — 9K and Ub = 3000 K. Stars show the zero approxi- 
mation -E(k) = h 2 eo(k)/2m. For the comparison with exper- 
iments, the dispersion curves are drawn to k max of the bcc 
lattice. 



It is the well-known Bogolyubov formula with the addi- 
tional multiplier 1/8 arising due to the boundaries (see 
I). Indeed, for each of the directions (1,0,0), (1,1,0), and 
(1,1,1), we need to continue curve (l52t to k max for a 
given direction. We obtained no anisotropy, but the gen- 
eral shape of the curve is proper: the phonon curve has 
no minimum for direction (1,0,0), a slightly pronounced 
minimum for (1,1,0), and a clear minimum for (1,1,1). 
This corresponds to experiments^ for the bcc lattice of 
He . The minimum arise due to a displacement of k max 
to the side of large k and to the presence of a minimum 
on the liquid-like curve E(k) — h 2 eo(k) /2m. To describe 
the curve more exactly and to obtain the anisotropy, we 
need to consider the next corrections to the equations. 

It is of interest that the simple formula (l52t describes 
the curves qualitatively correctly. This formula is uni- 
versal and is valid for a gas, a liquid, and a crystal, since 
these three states of substance are described by WFs with 
the same structure ©, (|57|) (the difference is only 
in the values of ki). 

The theoretical dispersion curves for a crystal and He 
II (see I) correspond to the experiment at Ub ~ 10 3 K. 
For Eq, the agreement with experiment for a liquid and 
a crystal holds at Ub ~ 10 2 K (see Sec. 2). However, it 
is assumed that the best potential for He 4 atoms is the 
Aziz potential 36 with much more larger value of £4, about 
2 x 10 6 K. This disagreement is related to the neglect of 
higher corrections or (more probably) to the fact that the 
potential at small distances has only an effective mean- 
ing. Therefore, the values of Ub can be very different for 
different processes. In some models, the agreement with 
experiment was attained with the Aziz potential. 



In Fig. [3] the dispersion curves are broken at k ma:c 
equal to a half of the minimum vector g of the recip- 
rocal lattice (for directions (1,0,0), (1,1,0), and (1,1,1) 
of a sc crystal, g = n(l/R, 0, 0), n(l/R, 1/R, 0) and 
it (1/R, 1/R, 1/R), whereas g is twice larger for the bcc 
lattice), since the dispersion curves for crystals are peri- 
odic with period g: 

£(k + g)=£(k). (53) 

The proof of this fact for one-particle WFs 3 .! can be eas- 
ily generalized to our case of TV-particle WFs. However, 
it is true for the cyclic BCs. For a three-dimensional 
crystal, the BCs are quite different and are close to zero 
ones. However, the majority of quasiparticles are local- 
ized wave packets (this is supported by the fact that the 
theory of transfer was constructed for wave packets and 
agrees with experiment). For such packets, the trans- 
lational invariance with the period of a lattice holds, 
and the conclusions of the theorem are valid. As for 
the standing waves |37|) - (|40|) . they are the sum of eight 
traveling waves, from which the wave packets can be con- 
structed. 

We note that the density of phonon states at the quan- 
tization of k by law ([5]) is the same as that for cyclic BCs. 

In the literatures^, the dispersion curves of crys- 
tals are usually calculated in the approximation of small 
oscillations ([1]). 



IV. THE CONDENSATE 

In Ref. and recently^, it was shown that the ideal 
crystal in the ground state contains no condensate atoms 
with k = 0. However, the condensate is possible in the 
presence of vacancies and other defect a 17 i 20 ~ 22 i 24 ' 26 . The 
condensate can be a key factor for the explanation of 
NCIM. But the experiment does not confirm 4 ^ a conden- 
sate of atoms with k = 0. 

It follows from vf^ © that the ideal sc crystal does 
not possess off-diagonal long-range order, but it have a 
condensate of atoms with k = ki. It is easy to verify 
if we switch-off the interatomic interaction: in this case, 
the exponents in ty © become equal to 1, and we have 
the product of sines with k = ki, i.e., all atoms are in the 
condensate with k = ki. If the interaction is switched- 
on, the condensate is exhausted. For a fluid, the conden- 
sate is also determined by a factor before the exponential 
function: under cyclic BCs, this factor (0 e lprj ' | p= o = 1) 

3 

generates the condensate of atoms with p = 0. For the 
fluid in a vessel, the factor is the product of sines with 
k = ki (see I), and the condensate is on the levels with 
k = ki,3ki,5ki,..4 3 .. 

Let us calculate the condensate for the ground state of 
a sc crystal. Under zero BCs, the condensate is deter- 
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mined by the formula 43 



iVk„ = 8 
N V 2 



J dridr 2 p(ri, r 2 ) sin (kn a xi) sin (fcn x x 2 ) 



x sin(fc„ ij 2/i)sin(fc„ B 7/2)sin(fc„ i zi)sin(fc„ i Z2), (54) 

p(r a ,r b ) = V J dr*^(r a ,f)*o(rb,r), (55) 

where n x ,n y ,n z are integers, n = (n x ,n y ,n z ), and r 
marks a collection of vectors r2, ...,rjv. Relation (|5J) 
yields 

p(r ,r 6 ) = sin (ki x x a ) sin (ki y y a ) sin (fcj^z,,) 
x sin (fc ix x&) sin (fej y y 6 ) sin (k lz z b ) 

x |A| 2 1/ J dr [exp{(5«(r a ,f))*} (56) 
x exp{s*W(r b ,r) + S* b (r a ,f) + 5 6 (r 6 ,f)} 

JV 

x J|{sin 2 (fe^Xj) sin 2 (fc^y, ) sin 2 (fo^)} 
i=2 



In ©, we consider only the one-particle part S\ 
Then 

p(r ,rfc) = sin^Xo) sin {h y y a ) sin (fcj^Za) 
x sin(fc; x Xb) sin (ki y b ) sin(ki sl z b ) 



(i) 



exp l ^p(qx) + (a: -> 2/, ^) >.(57) 



i? p (qj = -L ((5W(qJ)*e-^- + sf (qje^*) . 

(58) 

We use approximation (|28|) . (|26|) . Then is real, and 
we can replace e~ lq " x — > cosq^x and can sum only over 
> 0. We obtain 

p{r a ,r b ) = p^f x (x a )f x (x b )fy(y a )f y (y b )f z (z a )f z (z b ), 

(59) 



f x (x) = sin (ki m x) 



expj-^ £ ^(qDcos^r^j- 

(60) 

Since q x es = 2fc; x , 4fc; x , 6fc; x , . . ., we have 

f x (x) ^ sin(ki :c x)e {92COs{2k '- x)+9iCOs( - ik '- x ' )+ - } , (61) 

where g = 2S^ (jhJ x )/VN ~ 1. Function (J5TJ co- 
incides with ([2H1), d^fj). Substituting ([59]) and (jUJ in 
(1541 and expanding (|6T|) in a series, we obtain that the 
condensate levels with ~ A~ correspond to the wave 



vectors k = ((1 + 2j x )k lx ,{l + 2j y )k iy ,(l + 2j t )ki,) with 
j X ijyi3z — 0) 1) 2, 3, . . .; i.e., to the vector ki = k c and to 
larger vectors with odd multiple components. 

The distribution of atoms over levels depends on the 
barrier height U b , here we have the interesting pic- 
ture. The numerical calculation for He 4 atoms with 
R = 3.6 A and potential (03]) with a = 2 A, b = 4 A, 
U bd = -9K, and realistic values 500 if < U b < 5000 K 
gives 52 ~ lj 34 ~ —0.1, and the subsequent gj are 
small. With such gj, we obtain N^—^ ~ (1 — 92/2 + 
g' 2 /^) & p c 00 N/S. For the higher levels, we have: A^ k=3k sa 
Oft/2 - .g4/2) 6 pSo^/8, iVk= 5kc « (ff 2 2 /8 + 54 /2) 6 pL^/8, 
■^ r k=7k e ~ (525 l 4/4) 6 pJx 3 A^/8, etc. Hence, the levels 
with k = k c ,3k c and the intermediate ones with k = 
(k cx ,3k cy ,k cz ), (k cx ,3k cy ,3k cz ) (and with permutations) 
are filled. The rest levels are almost empty. Thus, we 
have 8 levels with N k ~ (g 2 /2) e p c 0O N/8, ~ p^N/8 3 atoms 
on each of them, in the sum N c ~ p^N/64. For com- 
parison, A^ kl ~ paoN/2 for a liquid^. It is easy to verify 
that 



Li; 



IP 3 



p 



dx 

L x 



[/x(x)] 2 



(62) 



With the obtained gj, we have P ~ (1— <72 + <?!)/2 — 1/2, 
i.e., p c x ~ 8/5^. The value of is given by the for- 
mula for a fluid. Let us consider the "crystal" cor- 
rection A 2 (k) (J20J) in Eq. dHJ for o 2 (k). Then, for 
C7& ~ 10 2 -10 4 K, the values of a 2 (k) decrease twice (in 
modulus) on the average in the significant interval of k 
(where a 2 is not small). In the zero approximation 44 , 



We know from experiment 



/ r 1 r-i a 2(k) 1 

Poo — ex P [-7? 2^ l-2a 2 (k)J 
k^O 

that Poo ~ 0.07 for He II. At twice less 02 (k), we obtain 
for He 4 atoms p 1 ^ ~ 1/3 and p c ~ 8/3. Hence, the 
condensate levels of a sc solid He 4 contain N c ~ 0.04A?~ 
atoms. This is only a rough estimate. 

At U b > 10 4 -10 5 K, we have 52 ~ 1, 34 takes val- 
ues from -2 down to -10, gj>6 ~ 0; the levels with 
k = 3k c ,5k c ,7k c and the intermediate ones are filled. 
At Ub > 10 6 K, the values of ge, gs, and 310 be- 
come large, which generates the condensates with k = 
9k c , like, 13k c . 

As is seen, the condensate structure in a crystal is sim- 
ilar to that of a fluid in a vessel 4 ^, but with the vector 
k c instead of ki and with a different distribution over 
levels. Within our method, it is impossible to calculate 
the condensate with high accuracy, because the signifi- 
cant corrections were omitted in almost all equations. In 
this case, the Monte-Carlo method can be efficient^. 

The structure of ([2]) helps us to imagine the conden- 
sate: ^0 contains N identical standing waves (the prod- 
uct of sines) resting on the walls by their wings. These 
waves form a single resonance classical wave, which mod- 
ulates the motion of atoms. Therefore, many atoms are 
characterized by the wave vector ki of this wave. It will 
be discussed in Sec. 8 that this wave can be considered 
as a particular kind of longitudinal sound. Therefore, the 
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condensate of N "zero-phonons" is present in the ground 
state. Thus, the resonance wave forms a crystal and sup- 
ports the condensate of atoms. In a fluid, the condensate 
is also formed by a wave, but with k = ki (see I and 
Ref.lH). 

Wave function ([2]) describes a simple rectangular lat- 
tice. The more complicated lattices can be constructed, 
by setting a relevant bare WF \fr bare instead of the prod- 
uct of sines \E r *£ re ®. The important point is whether 
\ff bare reflects the structure of a Wigner-Seitz cell, i.e., 
whether it equals zero on its surface. This is true for 
the sc lattice, but is not necessarily for other lattices. If 
ybare ^ q on the cell surface (we denote ^ b c % e ^ 0), then 
it is easy to guess the form of \fr bare for the bec and fee 
lattices: 

N 

*t = t"W2) J] { S m(k lv ( Xj - ai /2)) 

j=l+N/2 

x sm(ki y (yj - aj/2)) sin (fcj, (Zj - a;/2))} 

+ permutations. (63) 

*/cc e = ^7 e {N/A)^ x z v {NlA)% x {N/2)^ v x z (3A/4) 
+ permutations. (64) 

No+N/4 

K V (N Q ) = H {sm(v^*i.i J )sm(V2M i ) 

j=N + l 

x sm(ki z [zj +ai/2])}, (65) 
x.j = (xj — ai/2) cos (tt/4) + yj sin (7r/4), (66) 



It is of importance to confirm the existence of the con- 
densate in experiments. A neutron or a photon are scat- 
tered with the creation of a quasiparticle or elastically. In 
the second case, the wave vector of the scattering atom 
either is not changed or changes only the direction (to the 
opposite one), and the state of the crystal is invariable 
(except for the recoil). If the condensate is present, then 
the scattering with the momentum transference 2k c (or 
with a change in one-two components of the vector 2k c ) 
gives the intense peak. At = and at ^Itu ¥= 0, 

these changes in the momentum are equal to the vector g 
of the reciprocal lattice, and the scattering corresponds 
to the Bragg-Wolf (BW) peaks. But if 4-^ e = 0, then 
the number of condensate peaks is less than that of BW 
peaks. At $<!£jj e 7^ 0, the intensity of all BW peaks must 
increase at small T (in the presence of a condensate) pro- 
portionally to the amount of the condensate. For bare 
WFs ©, ([S5 ]l -(|55 ]) . the BW peaks (noncondensate ones) 
can be, apparently, interpreted as a result of the scatter- 
ing by zero-phonons with k = k c . 

For the bee, fee, and hep lattices, the Wigner-Seitz 
cell is complicated. Moreover, at ^eif = 0' ^ s crc " 
ation requires a wave with complicated structure. But, 
at ^ 0, the lattice can be formed from waves with a 

simpler structure of the type (|6"3"1) - (|6"6")) . Since the Nature 
uses simple structures, we suppose that 'I'cezT 7^ 0- 

We note that k = k c + <5k for the higher condensates, 
and the components 5k are multiple to the components 
2k c . Therefore, all condensates must be revealed in ex- 
periments as a part of the "base" condensate with k = k c . 



V. COMPARISON WITH THE TRADITIONAL 
APPROACH 



Vj = -(xj - ai/2)sin(7r/4) + yj cos (tt/4). (67) 

Here, the permutations symmetrize the WF to the Bose 
form, ai is the period of a lattice, and four functions in 
(|64| on the right-hand side depend on the coordinates of 
atoms with the numbers 1, . . . , A/4; A/4 + 1, . . . , A/2; 
N/2 + 1, . . . , 3N/A and 3N/4 + 1, . . . , N. These WFs 
imply that, for the bee lattice, k c = -K(i x + i y + i z ) / 'a\, and 
the condensate k are the same as that for the sc lattice. 
For the fee lattice, there are four composite condensates: 
with k c = 7r(i x + i y + i z )/a t and k^, = k(\[2\± + \[2\ y + 
i z )/ai with permutations of x, y, z. 

At **e« e = °' k c are different. For the sc lattice, the 
basis vectors of the reciprocal lattice are bi = 2wi x /ai, 
b 2 = 2iri y /ai, and b3 = 27ri z /et/, and the relation k c = 
So D /2 = (bi + b 2 + b 3 )/2 holds. For the bee lattice, 
we have bi = 27r(i y + i z )/ai, b 2 = 2ir(\ x + i z )/ai, and 
b,3 = 2tt(i x + i y )/ai, and the relation k c = gf® /2 — 
2tt(i x + i y + i z )/ai. 

We note that the solution for a one-particle WF with 
spherical orbit was studied in Ref. [23. There, it was 
assumed that the atoms of the condensate have k ~ 
2.1 A . This is close to our results. 



Let us compare the wave solution ([2J with the tradi- 
tional one (fTJ. In the last case, it is assumed that the 
atoms carry on small oscillations near lattice points. 

The traditional approach involves several assumptions. 
1) The lattice is set "by hands", it is not obtained from 
the Schrodinger equation. 2) It is assumed that the prob- 
ability density maximum in a cell is located at the lattice 
point, near which the atom carries on random oscilla- 
tions (atom is fixed likewise on rubber string, which role 
is played by the function ip (r — R)). Here, one more 
courageous assumption is hidden: that the mechanism of 
appearance of a lattice does not influence the motion of 
atoms in it. 3) It is assumed also that the boundaries 
have no effect on the solutions. Therefore, the realistic 
BCs close to zero ones are replaced by the cyclic condi- 
tions. 

Solution @ has no above-mentioned drawbacks. It is 
found with the use of natural zero BCs. The lattice and 
the probability distribution in a cell do not postulated, 
but they follow from the solution of the Schrodinger equa- 
tion. It is seen from the probability distribution that the 
motion of atoms is strongly affected by the mechanism 
(wave one) of formation of a crystal. In addition, the 
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observable quantities, in particular Eq and the energy of 
phonons, depend significantly on the BCs. This is caused 
by the fact that the walls change the Fourier expansion of 
the two-particle potential (see I: formulas (15)-(17) and 
Sec. 5) and affect the phonon frequency through ^o- Vi- 
sually, it is because the excitations of the quantum sys- 
tem are a standing waves rather then a particles; a wave 
keeps the memory about the wall and its length is mod- 
ulated by the wall. Such influence exists only for natural 
long-range potentials. In the above equations for E and 
E(k), this influence is manifested in the factor 1/8 at the 
potential Vs{k) and in the summation over the wave vec- 
tors multiple to tt/L, rather than 2n/L, in a corrections. 
To the effect of boundaries, we can refer also the correc- 
tions Qj (changing the phonon frequency E(k)/K), which 
arise from the product of sines in ^ and are induced by 
the interaction of a phonon with N zero-phonons of the 
ground state. 

Relation @ yields easily the traditional solution ([1}. 
Let each atom be localized near a lattice point, one atom 
per point. The points Rj correspond to the extrema of 
a sine. By expanding the sines in a Fourier series near 
points, the terms linear in Vj — Rj are absent, and the 
small quadratic terms can be taken up in the exponent. 
Thus, for the sc lattice, we have 

N 

]^[{sin (h x Xj) sin (ki y yj) sin (h z Zj)} 

3=1 

N 

-> ±JJ e -« 3 (ri-Ri) 2 /2 ) (68 ) 

3=1 

where a = ki T = ir/R (the close estimate a ~ 1 — 2 A 
was obtained in Ref. [j] from the other reasoning). By 

setting S$ = in ©, we reduce © to ©. 

Thus, the traditional solution ([I]) is a simplification of 
the wave solution @. It is significant that, in this case, 
the wave character of the solution is lost; but since the 
counting-off is made from the equilibrium positions of 
atoms, it is sufficiently simple to calculate 1 ^ -32 ' 35 ' 3 ^ 39 
(with several fitting parameters) Eq and E(k). In this 
respect, approach ([T]) has certain advantages, since the 
wave approach (|2|) yields a chain of complicated equa- 
tions with many large corrections for Eq and E(k). These 
equations cannot be solved exactly. 

The frequency of phonons E(k)/H is usually calculated 
in the harmonic approximation, where a phonon is a wave 
arising at a small deviation of atoms from equilibrium 
positions. However, it was noted above that the phonon 
frequency is affected by zero-phonons. In the language 
of oscillating atoms, this means that the atoms in the 
ground state are not in rest at lattice points (as is consid- 
ered in the harmonic approximation), but they oscillate 
intensively due to the motion in the field of zero-phonons. 
In this case, the potential energy is minimum, proba- 
bly, at lattice points. But it would be wrong to identify 
the deviations from "equilibrium positions" with the co- 
ordinates of real atoms, since, in this case, the intense 



zero oscillations would be lost. Therefore, the harmonic 
approximation is applicable only for the description of 
long-wave oscillations, when a crystal can be considered 
as a continuum. This approximation cannot be used for 
short-wave phonons, and the corresponding approaches 
in solid-state physics should be reconsidered, in our opin- 
ion. The agreement with experiment of phonon disper- 
sion curves calculated in the harmonic approximation in 
a number of works seems to be accident or is due to the 
choice of parameters. 

WF ([1]) is not a solution of the Schrodinger equation. 
The attempt to determine ip(r — R) in ([T]) from the 
Schrodinger equation leads 1 - to the solution ip(r — R) ~ 
e lkr , which is nonlocalized and, therefore, unphysical. It 
was assumed 1 - that this liquid-like solution arises due to 
the truncation of the cluster expansion. However, instead 
of e* kr , we can take a linear combination of exponents 
ip(r — R) ~ sin (ki^x) sin (ki y) sin (ki z z). The we arrive 
at the wave solution @, which is quite physical and sets 
a lattice with localized distribution of atoms. 

Recently, N. Prokof 'e^ considered a wave solution of 
the form 

N ( N L \ 

* = n(^X>( R ;- r *))> ( 69 ) 

Nl 

which contains a condensate in the state ^ ^(Rj — 

L 3=1 

r). However, this solution was recognized unphysical 3 , 
because the number Nl of the lattice points can differ 
from the number TV of atoms. Hence, the solution for a 
crystal is only one of the huge number of solutions, and 
the probability of its realization is too small. We note 
that N L in {2} can also be different from N. However, it 
is easy to see (see the following section) that we have no 
problems in this case. 

Thus, solution (0) can be obtained long ago. 

In the last time, the following new approaches to the 
description of quantum crystals are developed: Path In- 
tegral Monte-Carlo method 41 , variational Shadow WFs 
method 45 using the bare WF ^ and many fitting pa- 
rameters, and Shadow Path Integral Ground State pro- 
jector method 2 ^. The last method used the Bijl-Jastrow 
function, and the lattice arises due to the spontaneous 
breaking of symmetry. 

VI. GROUND STATE — LIQUID OR CRYSTAL? 

It is accepted that most substances in the ground state 
are crystals. Apparently, the ground state corresponds 
always to a liquid. 

The question about the structure of the ground state 
can be answered, by seeking a minimum of the energy. 
Consider a sc crystal. Its ground-state energy is given 
by formula (|36[) . In (|36[) . we fix the number of atoms N 
(i.e., the value of R) and change the number of points Nl 
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(i.e., the value of k c ). In the zero approximation (|32j) for 
a 2 (k), the curve E (k c ) is a parabola (see Fig. 4), where 

the point k c pa 1.41 A corresponds to JVx, = JV. If 
we descend somewhat downward along the curve, we ob- 
tain the states with Nl < N, where some cells have two 
atoms. In such cells, the distance between two atoms 
must be small, ~ A. Since the atoms have the almost 
hard core with a radius of 1.2 — 1.3 A, they strongly repel 
each other at distances of < 2.5 A. Therefore, such solu- 
tion is unstable: one of the atoms leaves the cell and will 
walk in the crystal, until it approaches the crystal sur- 
face. There, the atom evaporates or becomes fixed. As 
a result, the system transits in the state with Nl = N. 
If we lift upward along the curve, we obtain the states 
with vacancies (Nl > N). Such states have higher ener- 
gies. Hence, the stable state of the crystal with minimum 
energy has Nl = N, and each cell contains exactly one 
atom. This reasoning answers the objection in Ref. 0. 




Fig. 4: Function Eo(k c ) for a sc lattice of He 4 atoms at the 
fixed R = 3.85 A (rhombs) and the variable R under the con- 
dition Nl = N (circles). Stars show the underliquid at var- 
ious temperatures. We use the potential (|33[l with a — 2 A, 
b = 4k, U bd = -9K and U b = 81.2 K. 

The curve of rhombs indicates that E$(k c = 
1.6 A" 1 ) pa E {k c = 1.41 A" 1 ) + 1.1 K. At k c Pa 1.6 A" 1 , 
the system has 0.46 vacancy per atom. Hence, the vacan- 
cion energy E v pa 2.4 K, which is comparable, by the or- 
der of magnitude, with experimental value-- E v ~ 10 K. 

The dependence Eq(N) (circles in Fig. 4) is quite rea- 
sonable: the most favorable state of the lattice corre- 
sponds to the minimum (k c pa 1.41 A , R sa 3.85 A). 
By passing to a rarefied system (i? — > oo,k c — > 0), we 
have Eq — > 0. 

The minimum of the parabola of rhombs corresponds 
to the least possible fc c = fci = y/3n/L, at which WF @ 
transits in *o of a fluid, and E ([36]) is determined only 
by the term Eq for fluids. Hence, the ground state of a 
fluid (E w — 8.6K) lies below the ground state of the 
crystal (Eq pa — 1.27K). The states on the parabola near 



the minimum (e.g., with k c ~ 0.1 A ) correspond to a 
fluid partitioned into domains with many atoms. Such 
states must be unstable. 

We have studied Eq. (1551) obtained in a rather crude 
approximation. Its solution can be changed with regard 
for corrections and the exact potential (with a higher 
barrier). It can turn out that the ground state includes 
the admixture of vacancies^, though this is improbable. 

The main question is as follows: Is Eq of a crystal 
greater than Eq of a liquid or less, for the exact solu- 
tion? The answer is hinted by the node theorems!: "if 
the eigenfunctions of a self-adjoint second-order differen- 
tial equation are arranged in some region G under any ho- 
mogeneous BCs in the increasing order of the appropriate 
eigenvalues, then the nodal manifolds of the n-th eigen- 
function ty n divide the region G into at most n subregions 
for any number of independent variables" . The theorem 
was proved for WFs of the general form ^(ri, . . . , tn). 
We deal with the WFs of N identical Bose particles, and 
the conditions of the theorem are satisfied. We know that 
WF ([2]) with ki — > ki corresponds to the ground state 
of a liquid and has no nodes. By the theorem, namely 
this function describes the ground state of a system of 
Bose particles in a box. Hence, all remaining states, in- 
cluding the ground state of the crystal, have nodes and 
correspond to higher energies. But it can be so that the 
equations have no solutions for ^S/q corresponding to a 
liquid. Does of the crystal without nodes exist? It is 
clear that, at the wave structure of ([2]), the WF of the 
crystal has necessarily nodes. Can a nonwave structure 
be realized? The WF must turn to zero at the surface 
of the crystal and be symmetric under permutations of 
atoms. It is easy to see that the nonwave solution (fT]) 
is impossible without nodes: it does not feel the bound- 
aries and can become zero only at infinity. We cannot 
strictly prove this assertion, but we sure that only the 
wave solution satisfies the zero BCs. This implies that 
1) the ground-state energy of the crystal is always higher 
than that of a liquid composed of the same atoms; 2) the 
solution for describing a liquid exists necessarily for 
Bose atoms of all sorts (if the solution is absent at some 
density p of a liquid, it is possible to enter a region, where 
a solution exists, by varying p). 

However, the majority of substances in the Nature be- 
come crystals at low temperatures. The transition in the 
crystal state allows the system to decrease the energy 
by jump and is favorable. But it is seen in Fig. 4 that 
the principal minimum corresponds to a liquid. Appar- 
ently, by "having fallen" in the crystal state, the system 
cannot already pass to a deeper minimum corresponding 
to a liquid. To make this, the system must overcome a 
band of unstable state (or the energy barrier that can 
appear at the exact solution). Therefore, the majority of 
substances at low temperatures are crystals. 

The state of a liquid in the deep minimum can be called 
the underliquid (UL, stars in Fig. 4). Most probably, the 
majority of ULs has the superfluid phase. The tempera- 
ture (determined by quasiparticles) of UL can be equal to 
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the temperature of a crystal, but the total energy is less 
than that of the crystal. He 4 atoms have a large ampli- 
tude of zero oscillations. Therefore, the lattice is locally 
unstable, apparently, and transits in the UL state. In 
other words, He II is the single example of UL among 
inert elements. Possibly, other substances can be also 
transferred in the UL state (see Sec. 8). The amorphous 
bodies with microstructure of a liquid are not related, 
apparently, to UL, since the amorphous state is caused 
by a strong anisotropy of molecules, whereas we have 
considered the systems of spherical molecules. 

We arrive at a significant conclusion that a finite sys- 
tem of Bose particles of any sort (He, Ar, Ne, etc.) in 
the lowest state is a liquid, rather than a crystal, as is 
commonly accepted. Sometimes, the third law of ther- 
modynamics and the entropy-based arguments in favor 
of a crystal are discussed. However, the entropies of a 
crystal and a liquid in the ground state are identical and 
are equal to zero: S = ks In N s = ks In 1 = 0. In other 
words, the degrees of order of a fluid and a crystal are 
identical in the ground state, though a crystal seems vi- 
sually to be more ordered. 



VII. NATURE OF THE SUPERSOLID PHASE 

Almost all researchers arrived at the 
consensus2iS2r^2jSSiiiiiS of that the supcrsolid phase 
and NCIM 4 -^ are related to defects of the lattice. How- 
ever, the attempts to identify a carrier of NCIM with a 
specific defect met difficulties, which is not surprising. 
How can a crystal contain so many defects (or so 
extended defects) at ultralow temperatures T ~ 0.02 K 
that they connect 20% of atoms of the lattice, by 
ensuring the experimental value p s ~ 0.2/3? In our 
opinion, it is improbable. We note that the defects are 
an analog of quasiparticles; for comparison: the amount 
of quasiparticles in He II at T = 0.02 K is so small that 
they provide p n ~ 1O~ 7 ,0. Therefore, we suppose that 
the carriers of the effect are atoms of the ideal lattice. 

Let us consider various possibilities in detail. We start 
from vacancions. Since any crystal is in the gravity field, 
the vacancions undergo the action of the force directed 
upward. If the gas of vacancions is superfluid, then the 
vacancions should float up rapidly and evaporate from 
the surface. Under torsional oscillations, they must be 
transported to an internal surface of a crystal. Since this 
reasoning is valid for any massive defects, NCIM is not 
related to vacancions. 

The models of dislocation glass^,, superglass 2 - 1 , and 
grain boundaries 22 assume the existence of the conden- 
sate of atoms with k = 0, but it was not observed 42 : 
no < 0.003. In the model of screw dislocation network 24 , 
the superfluid component consists of atoms of the nuclei 
of dislocations, whose amount is obviously much less than 
10%. Therefore, this model does not explain the observed 
large values of p s ~ (0.1 — Q.2)p. The models of dislo- 
cation network met the analogous difficult y 18 ' 23 . As for 



the mechanism of grain boundaries 22 , it does not agree 
with the fact that NCIM was observed in a monocrystaliS. 
(where any grains are absent). 

As carriers can be atoms of the lattice. It was pro- 
posed already 2 ^, but no mechanism of correlations (of 
the condensate) was indicated. We propose the following 
scenario. For the sc lattice, = on the boundary 
of a Wigner-Seitz cell. Therefore, the atom cannot pass 
from a cell to another one. However, for the bec, fee, and 
hep lattices, $0 most probably on the boundary of a 
cell (see Sec. 4), and the tunneling of an atom from cell 
to cell is possible. Consider a bee crystal. According to 
Sec. 2, the atom is located near the orbit, whose size is 
twice less than that of a cell of the crystal. The orbits 
of the atom at the center of a cubic cell and the atom of 
one of eight vertices of a cube touch each other. There- 
fore, the WFs of such atoms are strongly overlapped. If 
the condensate includes more than a quarter of atoms, 
then each atom is contiguous to two or more atoms of 
the condensate. In this case, many neighboring atoms 
of the condensate can be joined by lines. Along such 
closed lines, the atoms can flow by means of the simulta- 
neous tunneling (according to the structure of WF (|63p , 
one atom cannot flow through a crystal, since the ze- 
ros of a sine create the impermeable planes in a certain 
distance). Due to the indistinguishability of atoms, the 
atoms belonging to the condensate and the overconden- 
sate change places. Therefore, all atoms participate in 
the flow. Since the atoms of the condensate are corre- 
lated, they can move as the whole, by representing the 
flowing component of a crystal. In order that this compo- 
nent be superfluid, its excitations must satisfy the Lan- 
dau criterion. The dispersion curve of such excitations 
was apparently observe d 13 ' 14 , and it satisfies the Landau 
criterion. We assume that the superfluidity is possible 
only in the case where the concentration n c of the con- 
densate is higher than the threshold one (n c J = n/4 for 
the bee lattice). Sec. 4 indicates that such n c is possible. 

Let us turn to the optic-like dispersion curve s 13 ' 14 . 
They are very close, and we can have no doubts that they 
represent the same mode. In Ref. this mode was con- 
sidered vacancional. But the experiment 1 ^ with a poly- 
crystal showed that it disappears at 0.2 K < T < 0.6 K. 
Though, the ordinary vacancional mode must be present 
at all T; it shifts only, as the density varies 4 ^. In Ref. IbH . 
the mode is referred to the superfluid component, and we 
agree with it. However, it is the mode related to a crystal 
(rather than to a liquid^), because the minimum is just 
at the Brillouin zone boundary. According to our ap- 
proach, this superfluid mode must appear at the cooling 
of a crystal down to the temperature, where the tunneling 
of the condensate starts. It is Tncim, the temperature, 
where the effect of NCIM arises. This is in agreement 
with datai^. In Ref. [H, the mode was observed also at 
T = 0.6 K, which is more than the known Tncim- But, 
in that case, a monocrystal was used. For it, Tncim is 
unknown and can exceed 0.6 K. We mention one more 
optic-like mode^S observed for the bee lattice of He 4 . It 
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should be noted that the data in Refs.tl3jfl4j|49l say noth- 
ing about the nature of the superfluid component, which 
can be vacancional. 

We note that since the flow can occur only through lat- 
tice points, and some condensate chains can break, the 
atoms of the condensate must be involved in torsional os- 
cillations of a crystal only partially. So that NCIM is pro- 
portional to qn c , where q < 1. The value of q is decreased 
by various defects. Therefore, q must strongly depend on 
the conditions of an experiment (this is observed), and 
p s can be small at high n c . The maximum of q for the 
crystals with a large ratio of area to volume^ can be ex- 
plained by the disappearance of large defects hampering 
the flow of the condensate. It is clear that the effect is 
maximum for the ideal crystal without admixtures and 
defects. 

We note that the destruction temperature for the con- 
densate must be significantly higher than TVc/m, since 
the tunneling flow requires n c > n^. r . 

The proposed tunneling mechanism agrees with the 
absence of the superfluidity of crystals at a pressure 
gradient^. According to Refs. [7||8|, the superfluidity 
should be accompanied by a deformation of the lattice 
(if defects are absent). But, at the tunneling mechanism, 
the atoms can move in the lattice without any change 
of its shape. Namely the torsional oscillations^ present 
the conditions for such motion. At the tunneling mech- 
anism, the crystal cannot also shift as a whole, since, to 
realize it, the surface atoms must tunnel outside of the 
crystal, where there are no lattice points. For the same 
reason, one more possible mechanism^ such as the leak- 
age of atoms of a fluid through a crystal is forbidden. All 
this corresponds to the conclusion^ that the mechanisms 
of flows for a crystal and a fluid are quite different. 

A decrease of NCIM after the annealing 6 can be related 
not to a decrease in the number of defects, but to the ap- 
pearance of lengthy defects. It is well known that a bullet 
piercing glass makes only a small hole, whereas a stone 
creates additionally a network of long cracks. In other 
words, the slow processes are accompanied by extended 
deformations of a crystal. Under the annealing, a crystal 
is firstly strongly heated and then is slowly cooled. In 
this case, the majority of defects disappear, but the re- 
maining dislocations are ordered into a three-dimensional 
network^, which hampers the flow of the condensate to 
a higher degree than many disordered dislocations. 

The effect of NCIM was mainly observed in polycrys- 
tals. The tunneling of atoms in them is possible, since 
there are no slits between microcrystals. 

The experimentally foundii law Ap ~ T 2 (it is equiv- 
alent to C v ~ T) can be because the condensate density 
is close to the threshold value: n c ~ rfj . In this case, 
the dimension of the network of lines, by which the con- 
densate flows, is close to 1. 

The jump of the rigidity of a crystal*^ at T « Tncim 
means that the same factor affects the rigidity and 
NCIM. We assume that this factor is the appearance of 
the superfluid component. As a result, a part of energy 



is transferred to modes of the superfluid component, and 
there occurs a redistribution of oscillatory modes of the 
system. This causes the jump of the rigidity. An in- 
crease of NCIM££ due to He 3 atoms can be related to 
the fact that He 3 atoms join dislocations and, under the 
action of the inertial force at torsional oscillations, turn 
a dislocation or move it to the surface of a crystal. 

It is necessary to understand why the heat capacity 
peak is less for purer crystals, and the temperature Tc v 
of the peak is less than Tnci m and is independent of the 
concentration of the He 3 admixture*- (though Tncim 
depends strongly on it). The first property is similar 
to the annealing effect. We assume that the reason is 
the same: purer (on the average) crystals contain more 
lengthy defects. Consider now the second and third prop- 
erties. The "excess" of C v is related to the superfluid sub- 
system and increases approximately linearly^ at low T. 
As T increases, a part of chains, on which the superfluid 
flow is realized, is broken due to the approach of n c to the 
threshold and the freezing-out of defects. Let a half of 
chains be broken. If many alternative chains remain, then 
NCIM is almost not changed, but the number of modes of 
the superfluid subsystem (i.e., C v also) decreases twice. 
NCIM decreases sharply only if the number of chains be- 
comes so little that the atoms lose the ability to flow. 
Apparently, as T decreases, the superfluidity arises grad- 
ually: first, a small number of atoms of the condensate 
can flow. In this case, NCIM and the peak of C v are not 
related to the phase transition, T^cim > Tq^, and the 
peak of C v is independent of the He 3 admixture and is 
caused by the dependence of the number of phonons in 
the condensate network of atoms on the properties of this 
network. 

As is seen, the proposed model can explain the super- 
solid phenomenon, but remains many questions unsolved. 
It should be clarified whether ^>q becomes zero on the 
surface of a Wigner-Seitz cell for the bec and hep lat- 
tices. It is also necessary to know which and how many 
defects are present in a crystal under various conditions, 
and how each sort of defects influences the condensate 
and its fluidity. The key moment for the model is the 
experimental discovery of a composite condensate with 
k = k c , 3k c , 5k c . 



VIII. DISCUSSIONS 

It is seen from formula @ that a crystal is formed by 
a standing wave in the probability field. A fluid has an 
analogous wave, but with k = ki = (tt/L Xi Tr/L y , ir/L z ). 
In I, it is shown that the wave in a fluid is similar to N 
standing sound waves with k = ki . The same arguments 
are true also for a crystal. Therefore, we can assert that 
a crystal in the ground state has N identical standing 
longitudinal acoustic phonons with k = k c . They are 
particular resonance zero-phonons. By the structure of 
corrections, they differ from ordinary phonons described 
by WF f57 |) -(|lD l) . It is significant that these standing 
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zero-phonons create a crystal lattice. In other words, the 
periodicity of a crystal is caused by the periodicity of a 
sound wave. 

The number of resonance phonons is equal to the num- 
ber of atoms. At such huge occupation number, these 
phonons can be considered as a single classical sound 
wave. But if a crystal is a wave, we can try to control its 
state with the help of sound and electromagnetic waves. 
In particular, it would be possible to create or to destroy 
crystals. 

Well-known are the legends concerning N. Tesla£2, who 
induced vibrations of building's walls with the help of a 
small mechanical oscillator with the resonance frequency 
in the ultrasound region. Apparently, N. Tesla excited 
the eigenmodes with k ([5]). The particular resonance 
frequency for crystals is the frequency of zero-phonons, 
which is equal to the difference of Eq/H of a crystal and 
Eq/H of a fluid. The wave with such a frequency forms a 
crystal itself. 

In Sec. 6, the state of underliquid is predicted. Pos- 
sibly, this state can be obtained by means of the wave 
action on a crystal or a rapidly cooled fluid (so that 
the fluid will avoid the state of crystal and become UL). 
Through a liquid, it is necessary to transmit monochro- 
matic sound or electromagnetic waves with a wavelength 
comparable (but not equal to) with the period of a lat- 
tice. It is possible to use an x-ray or gamma-laser with 
0.5 A < A < 10 A. But such lasers have not been created 
till now, to our knowledge. A crystal can be destroyed, 
possibly, by a wave with the frequency of a zero-phonon. 
All this remains else at the level of speculations, but it is 
interesting to study these quastions. 

The crystallization of a liquid can be related to reso- 
nance phenomena in the system of phonons. The centers 
of crystallization are usually considered in the language 
of interacting atoms. But, according to the wave solu- 
tion, these centers are, most likely, growing wave packets 
with k = k c . 

We assume that the wave principle is the general one 
for the formation of crystals. As an important task, we 
indicate the search for the solutions for the bcc, fee, and 
hep lattices, which are most spread in the Nature. 



The properties of quantum crystals are fine and arouse 
the feeling of admiration. In addition, it seems clearly 
that Lady Science moves away from Truth sometimes in 
her walks. But She goes not so quickly, as the Universe 
expands, and can return always. 

The present work is devoted to the memory of Petr 
Ivanovich Fomin. 



X. APPENDIX. CALCULATION OF Ci(j) 

Formula (TT21 arises from the expansion of f(x) (Tn| in a 
Fourier series as 5 — > 0. Instead of f(x), we can take any 
function, which satisfies the conditions for the Fourier 
expansion to exist and passes to cot(fc; x a;) as 5 — > 0. 
Since the input function can be expanded in a Fourier 
series at 5 ^ 0, the quantity Ci x (j x ) can be determined 
by the formulas of Fourier analysis: 



CiM = -CiA-J x ) = lim / dxf{x)e-^' L * 
<5^o L x J 



-5- f dxcot(k l:c x)e- l2 ^ x/L * 
Li x J 
o 

1/2 



2i / dx sin (2irj x x) cot (ttI x x). 



(70) 



Integral (IT0"|) is defined, despite the discontinuities of a 
cotangent, and can be calculated in the sense of the prin- 
cipal value of an improper integral^. By induction, re- 
lation dTOJ) yields Ci x (j x ) = -i at j x = l x ,2l x ,..., and 
C; x = i at j x = —l x , —2l Xl . . .. However, the proof of the 
third part of (TT3l) , namely the vanishing of C; x at the rest 
j x , is not a simple task. 

We use the following trick. In I, it was shown that 

cot(Wi) = -*E(j/b1)^ /L - (71) 



IX. CONCLUSIONS 

Our most important conclusion concerns the wave na- 
ture of Bose crystals. The wave properties yield the con- 
densate of atoms and the possible superfluidity of a crys- 
tal. If the crystals in the Nature are created by standing 
waves, it is quite beautiful. 

The complexity of properties of quantum crystals is 
related to the fact that they have five subsystems: atoms 
of the lattice, atoms of the condensate, quasiparticles 
in both these systems, and various defects. This makes 
it difficult clarifying the nature of the supersolid phase. 
Though He II has only two subsystems, the superfluid 
subsystem and that of quasiparticles, the nature of its 
superfluidity was understood in more than one decade. 



This equality should be understood in the same meaning 
as (1121) , which was obtained in I from the Fourier expan- 
sion of the smooth function cot {k\ x x + 5\ x ) as 5\ x — > 0. 
Formula (|TTj) follows also from the sum of two geometric 
progressions: 



cot{irx/L) = < —i 



p j(-5+i)2irx/L 



3 = 1,2, 



+ i ' 

j=-l,-2. 



j(8+i)2irx/L 



(72) 



<5^0 



where S > 0, x e [0, L]. Replacing x — > l x x in ([TTj) or 
(1721 . we obtain (TT5]). 
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